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Executive  Summary 

The  main  goal  of  this  research  was  to  provide  fundamental  data  on  the  dynamic  behavior  of 
laminar  flames  for  a  wide  range  of  conditions.  The  parameters  considered  included  the  fuel  type, 
reactant  composition,  flame  temperature,  and  combustion  mode.  The  derived  experimental  data 
constitute  a  basis  for  partially  validating  the  combustion  chemistry  of  a  large  number  of  fuels 
ranging  from  hydrogen  to  gaseous  and  liquid  hydrocarbons  and  alcohols.  Emphasis  was  given  to 
kinetics  regimes  that  have  not  been  probed  systematically  in  past  flame  studies  and/or  in  which 
current  knowledge  of  combustion  chemistry  appears  to  be  inadequate.  The  understanding  of  the 
oxidation  kinetics  of  such  fuels  is  of  immediate  interest  to  the  design  of  advanced  air-breathing 
propulsion  devices.  The  experiments  were  performed  in  the  counterflow  configuration,  and  the 
flow  velocities  were  determined  through  the  use  of  Digital  Particle  Image  Velocimetry.  Laminar 
flame  speeds  and  extinction  strain  rates  were  determined,  as  they  are  both  fundamental  flame 
properties  and  sensitive  to  chemical  kinetics.  Additionally,  ignition  states  in  both  premixed  and 
non-premixed  configurations  were  determined  for  a  wide  range  of  fuels.  The  experiments  were 
modeled  through  the  use  of  a  quasi  one-dimensional  code  and  included  detailed  descriptions  of 
chemical  kinetics,  molecular  transport,  and  thermal  radiation.  The  results  of  this  research 
convincingly  showed  that  current  state  of  the  art  kinetics  models  have  notable  limitations  in 
predicting  flame  properties  even  for  simple  molecules  such  as  H2  and  C1-C2  hydrocarbons,  with 
these  limitations  being  more  profound  for  the  more  complex  liquid  hydrocarbons.  It  also  was 
found  that  for  H2  flames,  uncertainties  in  diffusion  coefficient  calculations  could  have  a  greater 
effect  on  the  predicted  flame  properties  compared  to  kinetics,  a  result  that  further  complicates 
current  approaches  aiming  towards  kinetics  validation. 

This  research  was  performed  by  the  PI,  two  Research  Associates  (Dr.  Yufei  Dong  and  Dr. 
Gurhan  Andac),  and  two  graduate  students  (Jeff  Langille  and  Adam  Holley).  The  following 
j  oumal  papers  were  authored : 

1.  "Measurement  of  Laminar  Flame  Speeds  through  Digital  Particle  Image  Velocimetry: 
Mixtures  of  Methane  and  Ethane  with  Hydrogen,  Oxygen,  Nitrogen,  and  Helium,"  by  Y. 
Dong,  C.  M.  Vagelopoulos,  G.  Spedding,  and  F.  N.  Egolfopoulos,  Proceedings  of  the 
Combustion  Institute  29:  1419-1426  (2002). 

2.  "Non-Premixed  Ignition  of  H2  and  CO  by  Vitiated  Air  in  Counterflow  Configurations,"  by 
J.A.  Langille,  Y.  Dong,  M.G.  Andac,  F.N.  Egolfopoulos,  and  T.T.  Tsotsis  submitted  for 
consideration  at  the  30*  Symposium  (International)  on  Combustion. 

3.  "Premixed  Flame  Extinction  of  Practical  Liquid  Fuels:  Experiments  and  Simulations,"  by 
A.T.  Holley,  A.  Balteria,  Y.  Dong,  Y.  Fan,  M.G.  Andac,  and  F.N.  Egolfopoulos,  submitted 
for  consideration  at  the  30*  Symposium  (International)  on  Combustion. 
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Dong,  A.T.  Holley,  M.G.  Andac,  F.N.  Egolfopoulos,  S.G.  Davis,  P.  Middha,  and  H.  Wang, 
submitted  for  consideration  at  the  30*  Symposium  (International)  on  Combustion. 
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Technical  Summary 

1.0  Introduction 

The  accurate  knowledge  of  the  oxidation  kinetics  of  hydrocarbon  fuels  is  essential  for  the 
design  of  the  next  generation  of  air-breathing  engines  operating  at  speeds  and  altitudes  that  are 
noticeably  greater  compared  to  the  ones  currently  achieved.  State-of-the-art  knowledge  of 
hydrocarbon  oxidation  chemistry  is  limited  to  that  of  H2  and  CH4.  For  example,  the  GRI-3.0 
mechanism  [1]  closely  predicts  a  wide  range  of  flame  properties  for  H2  and  CH4  mixtures  with 
air,  especially  at  high  temperatures.  However,  even  H2  and  CH4  oxidation  kinetics  at  low 
temperatures  have  not  been  probed  systematically  in  previous  flame  studies.  The  oxidation 
kinetics  of  higher  carbon  hydrocarbons  are  subjected  to  significant  uncertainties  even  at  the  C2- 
level  [e.g.,  2].  Furthermore,  systematic  flame  studies  of  liquid  hydrocarbons  are  limited 
compared  to  the  gaseous  ones. 

This  research  was  both  experimental  and  numerical.  The  experimental  measurements 
focused  on  the  determination  of  global  flame  properties  of  relevance  to  the  phenomena  of 
ignition,  propagation,  and  extinction,  and  a  number  of  innovations  were  introduced.  For 
example,  while  propagation  and  extinction  typically  are  treated  as  high-temperature  phenomena 
for  hydrocarbons,  subsets  of  their  kinetics  also  were  probed  at  low  temperatures.  This  probing 
was  achieved  by  studying  ultra-lean  and/or  highly  diluted  hydrogen  flames.  Such  flames  can  be 
sustained  at  temperatures  as  low  as  900  K,  a  regime  of  interest  for  hydrocarbon  ignition. 
Additionally,  the  high-temperature  kinetics  were  validated  against  flame  extinction  data  obtained 
for  a  wide  range  of  fuels,  such  as  CO,  C1-C3  hydrocarbons,  and  liquid  fuels  (  Cy-Cg  hydrocarbons 
and  alcohols).  While  laminar  flame  speeds  cannot  always  be  determined  directly,  extinction 
strain  rates  (in  counterflow  configurations)  can,  and  a  methodology  has  been  introduced  to  do  so. 
It  has  been  shown  that  similar  kinetics  pathways  control  propagation  and  extinction  [2].  An 
alternative  technique  also  was  introduced  for  studying  flame  ignition  in  counterflow 
configurations  by  using  vitiated  air  that  is  produced  from  the  oxidation  of  ultra-lean  H2/air  flames 
as  the  ignition  source,  rather  heated  air.  The  ignition  behavior  of  several  fuels  thus  was  tested. 
The  final  innovation  included  constant  flame  temperature  experiments  for  a  given  fuel  and 
pressure  and  variable  fuel  concentration.  These  conditions  could  be  achieved  by  controlling  the 
ratio  of  O2  to  N2  in  the  oxidizer.  Thus,  any  kinetics  deficiencies  of  sufficiently  off-stoichiometric 
mixtures  would  not  be  masked  behind  the  uncertainties  associated  with  small  values  of  measured 
properties.  Such  is  the  case  when  the  flame  temperature  is  also  reduced  with  stoichiometry. 
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2.0  Experimental  and  Numerical  Methodologies 

The  experiments  were  performed  in  the  opposed-jet  configuration  and  through  the  use  of  a 
recently  developed  correlation  image  velocimetry  (CIV)  technique  [3],  a  variation  of  the  digital 
particle  image  velocimetry.  Laminar  flame  speeds,  S°,  were  determined  at  the  point  of  transition 
from  a  planar  to  a  Bunsen  flame  as  had  been  proposed  by  Vagelopoulos  and  Egolfopoulos  [4] 
(VE98).  The  use  of  CIV  allowed  the  determination  of  the  instantaneous  velocity  field  just  before 
the  transition  initiates.  Thus,  the  instantaneous  strain  rate,  K,  and  the  attendant  5"  could  be 
measured  directly.  The  extinction  strain  rate,  also  was  measured  directly  by  establishing 
flames  very  close  to  extinction  and  by  reducing  (for  fuel-lean  conditions)  the  fuel  flow  rate  [5,6] 
rather  than  increasing  the  strain  rate;  the  latter  procedure  also  requires  extrapolations.  The 
experiments  were  modeled  by  using  a  quasi  one-dimensional  code  and  detailed  description  of 
molecular  transport,  chemical  kinetics,  and  thermal  radiation  [7].  The  chemical  kinetics  schemes 
tested  included  (1)  the  GRI30  mechanism  [1]  for  C1-C3  hydrocarbons,  (2)  the  H2/O2  mechanism 
of  Mueller  et  al.  [8],  (3)  a  mechanism  developed  by  Wang  and  coworkers  [9]  for  Ci-Ce 
hydrocarbons,  (4)  three  mechanisms  for  methanol  and  ethanol  oxidation  [10-12],  and  (5)  two 
mechanisms  describing  Cv-Cg  hydrocarbon  oxidation  [13,14]. 

3.0  Results  and  Discussion 

3.1  Propagation  and  Extinction  of  Cj-Cj  Hydrocarbon/Air  Premixed  Flames 

Figure  1  compares  the  measured  and  those  determined  in  VE98  with  laser  Doppler 

velocimetry  (LDV)  for  CHVair  mixtures  at  various  equivalence  ratios  (<|)’s).  Ref  3  contains  the 
pertinent  methodology  and  results  described  in  this  section.  The  two  sets  of  data  are  in  close 
agreement,  which  is  an  independent  validation  of  the  accuracy  of  the  ClV-derived  data.  The 
numerically  determined  also  are  shown  in  Fig.  1,  and  close  agreement  can  be  seen  in 
general.  The  inset  depicts  the  dependence  of  the  predicted  Si  on  the  number  of  grid  points  used 
in  the  simulations.  There  is  a  systematic  over-prediction  of  the  experimental  5"’s  as  the  mixture 
becomes  progressively  leaner. 

The  discrepancies  between  experimental  and  predicted  Si  for  very  lean  mixtures  have  not 
attracted  proper  attention,  as  the  observed  differences  of  1-2  cm/s  typically  are  considered  to  lie 
within  experimental  uncertainties.  This  difference  occurs  because  as  (j)  decreases,  S^^  decreases 
as  a  result  of  the  flame  temperature  reduction.  Thus,  if  Si  is  to  be  used  for  validating  chemical 
kinetics  for  a  wide  range  of  (|)’s,  the  proper  way  to  do  it  is  by  keeping  the  flame  temperature 
constant  as  ^  changes  so  that  concentration  effects  are  isolated  from  those  of  temperature  [15]. 

Figure  2  depicts  experimental  and  predicted  Si's  for  CH4/O2/N2  mixtures  as  functions  of  <j), 
keeping  the  adiabatic  flame  temperature,  fj,  constant  by  modifying  the  relative  O2/N2  free- 
stream  concentrations;  was  determined  by  equilibrium  calculations.  For  all  three  f/s 
investigated  i.e.  1900,  2000,  and  2100  K,  Si  exhibits  a  non-monotonic  behavior.  For  ([)<  1 
mixtures,  Si  increases  as  <))  decreases.  Species  consumption  and  sensitivity  analyses  showed  that 
as  <|)  decreases  the  rate  of  the  main  branching  reaction  H  +  O2  ^  OH  +  O  (Rl)  is  enhanced  as  the 
O2  concentration  increases  for  the  same  .  For  (|)  >  1  mixtures.  Si  first  increases  and  then 
decreases  as  ^  increases.  Si  increases  because  the  O2  concentration  increases  relative  to  N2  in 
order  to  maintain  the  same  Tad-  Thus,  the  rate  of  Rl  is  enhanced.  Subsequently,  the  mixture 
reactivity  is  reduced  as  the  CH4  concentration  increases  to  the  point  that  it  effectively  scavenges 
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H  radicals  through  the  reaction  CH4  +  H  ^  CH3  +  H2  (R2)  competing  for  H  radicals  with  Rl, 
whose  rate  is  reduced.  As  a  result  the  burning  rate  is  also  reduced. 


Fig.  1 .  Comparisons  between  the  experimentally 
measured  5"  data  for  atmospheric  CIVair 
mixtures  with  those  of  VE98  and  predictions 
using  the  GRI30  mechanism. 


Fig.  2.  Variation  of  experimental  and 
predicted  5"  with  <j)  and  flame  temperature 
for  atmospheric  CH4/O2/N2  mixtures. 


Comparing  the  experimental  and  predicted  S^s,  there  is  a  close  agreement  for  near- 
stoichiometric  and  rich  flames.  While  for  lean  flames  the  agreement  is  very  close  for 
7’^_^  =  2100K,  for  the  lower  T^^’s  the  experiments  were  notably  over-predicted  by  the 
simulations.  The  observed  discrepancies  are  greater  than  the  experimental  uncertainty,  which 
had  been  determined  to  be  less  than  3%. 

Extinction  data  were  derived  for  fuel-lean  mixtures  of  air  with  CH4,  C2H6,  C2H4,  C2H2,  CsHg, 
and  C3H6  by  utilizing  the  single  flame  configuration  that  results  by  counterflowing  mixtures 
against  an  air  jet.  Given  that  it  was  possible  to  measure  the  values  directly,  focus  was  on 
using  for  model  validation.  Because,  there  is  always  a  small  but  finite  uncertainty 

associated  with  the  determination  of  ’s,  either  as  a  result  of  the  presence  of  residual  stretch  or 
the  need  of  using  fitting  formulas  to  extract  5“.  Figures  3  through  8  depict  experimental  and 
numerical  for  single,  lean  premixed  CHU/air,  C2H6/air,  C2H4/air,  C2H2/air,  C3H8/air,  and 

CgHe/air  flames.  The  inability  of  the  GRI30  mechanism  [1]  to  predict  other  than  CH4  and 
C2H6  flames  is  apparent  and  in  agreement  with  a  previous  study  [2].  The  Wang  mechanism  [9] 
was  tested  against  the  C2H2/air,  C3H8/air,  and  C3H6/air  data.  Although  the  simulations  under¬ 
predicted  somehow  the  experiments,  the  agreement  is,  in  general,  fairly  good. 
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Equivalence  ratio,  ^ 

^  Fig.  3.  for  single  CHVair  flames 


Fig.  4.  for  single  flames 


Fig.  6.  for  single  C2H2/air  flames 


Fig.  7.  for  single  CsHg/air  flames  Fig.  8.  for  single  CsHe/air  flames 
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3.2  Premixed  Flame  Extinction  of  Mixtures  of Liquid  Fuels  with  Air 

The  extinction  of  mixtures  of  practical  liquid  fuels  with  air  also  was  studied  in  the  single 
flame  configuration;  Ref  5  contains  the  pertinent  methodology  and  results  described  in  this 
section.  The  introduction  of  the  liquid  phase  into  the  air  under  room  temperature  was  achieved 
through  the  use  of  an  evaporator  and  a  high-precision  liquid  pump.  The  liquid  was  atomized  into 
a  stream  of  heated  air,  and  subsequently  the  temperature  of  the  mixture  was  reduced  to  ambient. 
The  liquid  fuel  existed  in  the  vapor  phase  as  long  as  its  partial  pressure  was  maintained  below  its 
vapor  pressure  at  the  prevailing  ambient  temperature.  The  liquid  fuels  that  were  tested  were 
CH3OH,  C2H5OH,  n-C7Hi6,  and  iso-CgHig.  The  A'^^/s  of  mixtures  of  those  fuels  with  air  at 
ambient  temperature  are  shown  in  Fig.  9. 


Fig.  9.  for  mixtures  of  air  with  liquid  fuels 


The  ATgjc/'s  of  the  alcohol  fuels  are  notably  higher  than  those  of  n-CyHn,  and  iso-CgHig,  while 
there  are  minor  differences  between  the  measured  of  CH3OH  and  C2H5OH  flames. 

Furthermore,  the  n-C7Hi6  flames  are  more  resistant  to  extinction  than  iso-CgHig  flames.  These 
observations  are  important  in  terms  of  combustion  stability. 

Figure  10  depicts  the  and  5"  values  for  methanoFair  mixtures  at  various  (l)’s;  S^'s  were 
taken  from  Ref  16.  The  data  are  slightly  under-predicted  by  the  Held  and  Dryer  [10] 
mechanism  (HD98)  for  ^  <  0.75  and  the  agreement  is  much  closer  for  (|)  >  0.75.  Comparing  the 
S°'s  reveals  that  the  HD98  mechanism  over-predicts  the  data  for  the  entire  range  of  <|)’s.  The 
Fischer  et  al.  [11]  mechanism  (FDCOO)  satisfactorily  predicts  the  S°'s  for  lean  (f)’s;  it  under¬ 
predicts  the  ’s  in  the  same  (|)-range  by  a  factor  of  nearly  3.  Figure  11  depicts  and  5" 
values  for  ethanoFair  mixtures  at  various  (])’$.  ^'"'s  were  taken  from  Ref  17.  Both  the  FDCOO  and 
the  Marinov  [12]  (MRN99)  mechanisms  under-predict  the  experimental  Kg^t  by  nearly  a  factor 
of  2.  The  S^'s  are  slightly  over-predicted  over  the  entire  range  of  (|)  by  FDCOO.  The  MRN99 
mechanism  slightly  over-predicts  ■S'^'s  on  the  lean  side  and  slightly  under-predicts  them  for  near- 
stoichiometric  and  rich  flames. 
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Figure  12  depicts  the  and  values  for  n-heptane/air  mixtures  at  various  (l)’s.  The 
experiments  were  modeled  using  the  Davis  &  Law  [13]  mechanism  (DL98).  While  Si's  are 
closely  predicted  on  the  fuel-lean  side  and  under-predicted  on  the  fuel-rich  side,  the  opposite  is 
true  for  More  specifically,  on  the  fuel-lean  side  the  experimental  are  under¬ 

predicted  by  nearly  a  factor  of  2,  while  on  the  fuel-rich  side  a  close  agreement  is  observed.  The 
predicted  Kg^f  peaks  at  about  (()  =  1 .22,  while  the  experimental  Kg^^f  at  about  (j)  =  1.08. 

Figure  13  depicts  the  Kg^^  and  SI  values  for  iso-octane/air  mixtures  at  various  (|)’s.  The 
experiments  were  modeled  using  both  the  DL98  and  Pitsch  et  al.  [14]  (PPS96)  mechanisms.  It  is 
apparent  that  both  mechanisms  fail  to  predict  the  experimental  's,  especially  on  the  fuel-lean 
side  where  the  experimental  data  are  under-predicted  by  factors  of  2  to  3.  Furthermore,  the 
experimental  Kg^t  peaks  around  (|)=  1.15,  while  the  DL98  predicts  the  peak  approximately 
around  ([)=  1.25,  and  PPS96  predicts  the  peak  around  ([)=  1.30.  On  the  fuel-rich  side  PPS96 
over-predicts  the  experimental  ^ext  's>  while  DL98  under-predicts  them.  The  observed 
disagreements  on  the  fuel-rich  side  are  less  that  the  ones  observed  on  the  fuel-lean  side.  Both 
mechanisms  result  in  rather  similar  SI  values  for  the  entire  (|)-domain.  They  both  predict  closely 
the  experimental  Si's  for  fuel-lean  mixtures,  while  they  under-predict  them  on  the  fuel-rich  side. 
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Equivalence  Ratio, 

Fig.  12.  and  S°  for  n-heptane/air  flames 


Equiv'alence  Ratio,  <|) 

Fig.  13.  Kgxt  and  for  iso-octane/air  flames 


Comparisons  between  experimental  data  and  numerical  predictions  revealed  that  for  the 
liquid  fuels  studied  in  this  investigation,  kinetics  mechanisms  that  closely  predict  laminar  flame 
speeds  fail  to  predict  the  attendant  extinction  strain  rates.  However,  this  inconsistency  was  not 
found  to  be  the  case  for  mixtures  of  air  with  methane,  a  fuel  whose  oxidation  kinetics  are  well 
characterized.  The  observed  discrepancies  could  be  attributed  to  the  fact  that  the  validation  of 
kinetics  mechanisms  of  the  liquid  fuels  is  not  as  comprehensive  as  that  for  methane.  As  a  result, 
the  ability  of  such  mechanisms  to  predict  flame  phenomena  even  at  relatively  similar  conditions, 
such  as  propagation  and  extinction,  may  be  compromised.  This  finding  is  of  interest  for  large- 
scale  simulations  of  practical  combustors. 


3.3  Chemical  Kinetics  and  Diffusion  Effects  in  Premixed  Extinction  of  H 2/air  Flames 

The  oxidation  kinetics  of  H2  were  tested  under  fuel-lean  conditions  for  which  Ihe  presence  of 
thermo-diffiisional  instabilities  does  not  allow  for  the  determination  of  laminar  flame  speeds; 
Refs.  6  and  18  contain  the  pertinent  methodology  and  results  described  in  this  section.  Under 
such  conditions  the  flame  temperatures  are  well  below  those  encountered  in  hydrocarbon  flames. 
Such  temperatures  are  relevant  to  the  process  of  ignition.  The  extinction  strain  rates,  K^f,  were 
determined  and  used  to  validate  the  kinetics,  as  it  is  known  that  both  extinction  and  propagation 
are  controlled  by  similar  chain  mechanisms  [2].  The  presence  of  positive  stretch  suppresses  the 
formation  of  cellular  structures.  Figure  14  depicts  comparisons  between  the  experimental  and 
predicted  using  the  GRI30  [1]  and  Mueller  et  al.  [8]  ("MKYD99'')  H2/O2  mechanisms. 

Both  experiments  and  simulations  were  performed  in  the  single-flame  configuration  in  which  a 
mixture  jet  counterflows  against  an  ambient  temperature  air-jet.  The  agreements  are  not  close, 
and  there  is  a  difference  in  the  "slope"  of  the  results.  While  the  predicted  K^/s  by  the  two 
mechanisms  "merge”  at  high  temperatures,  they  deviate  significantly  from  each  other  at  lower 
temperatures.  The  GRI30  mechanism  appears  to  under-predict  the  K^^'s  of  ultra-lean  H2/air 
flames  by  a  factor  of  two  (!)  at  low  temperatures.  At  high  temperatures,  both  mechanisms  appear 
to  over-predict  the  experimental  K^^f's  by  nearly  50%.  The  numerical  simulations  that  are 
shown  in  Fig.  14  were  performed  by  using  a  large  number  of  grid  points,  N.  It  was  found  that, 
by  varying  N  from  300  to  2000-3000,  could  increase  by  as  much  as  25%,  reaching 
asymptotically  the  values  that  are  reported  in  Fig.  14  as 


9 


Fig.  14.  for  single  Ha/air  flames 

Through  detailed  sensitivity  analysis,  it  was  also  shown  that  the  extinction  states  of  such 
mixtures  are  not  only  sensitive  to  the  kinetics  but  also  to  the  species  diffusivities.  Thus, 
validation  of  models  against  flame  data  of  such  highly  diffusive  mixtures  must  be  done  by 
considering  the  uncertainties  associated  with  both  kinetics  and  diffusion. 

In  view  of  this  consideration,  an  additional  contribution  of  this  study  was  the  advancement  of 
an  improved  treatment  of  the  mass  diffusion  coefficients  [6].  Figure  17  depicts  comparisons 
between  the  experimental  and  predicted  for  these  models  (GRBO,  "MKYD99",  and  a  trial 

model  "TRM04"  and  an  optimized  model  "OPM04"  [19]).  The  simulations  shown  in  both  Figs. 
15a  and  15b  included  the  use  of  the  updated  multi-component  diffusion  model.  Significant 
disagreements  exist  between  experimental  and  predicted  Kgxt%  and  there  is  a  difference  in  the 
“slope”  of  the  results.  Detailed  simulations  also  showed  that  the  predicted  extinction  strain  rates 
based  on  the  updated  formulation  of  dif&sion  coefficients  were  different  by  ~15-25%  from  those 
obtained  using  the  Sandia  transport  coefficient  compilation  (Fig.  16).  Similarly,  it  was  found  that 
by  using  mixture-averaged  versus  full  multi-component  formulation  could  affect  the  predicted 
extinction  strain  rates  by  as  much  as  25-40%.  From  a  kinetics  point  of  view  differences  of  the 
order  of  15-40%  in  predicting  global  flame  properties  are  rather  large  and  could  have  a  major 
impact  on  rate  constants. 
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Adiabatic  Flame  Temperature,  K 
1190  1311  1428  1539  1190  1311  1428  1539 


Equivalence  Ratio,  ^ 

Fig.  15.  Experimental  and  numerical  for 
single  ultra-lean  H2/air  flames.  The  effect  of 
reaction  model  is  shown. 


Adiabatic  Flame  Temperature,  K 
1190  1311  1428  1559  1190  1311  1428  1539 


Fig.  16  Experimental  and  numerical  for 
single  ultra-lean  H2/air  flames.  The  effect  of 
diffusion  model  is  shown. 


In  order  to  moderate  the  diffusion  effects  additional  experiments  were  performed,  by 
determining  Kg^/s  for  N2-diluted  near-stoichiometric  H2/air  mixtures.  The  results  are  shown  in 
Figs.  17  and  18  for  (|)  =  0.777  and  0.965,  respectively.  In  these  studies  the  associated  Tad's  vary 
from  1600-1700  K  down  to  ~1200  K,  thus  accessing  similar  flame  temperatures  as  those  of  ultra¬ 
lean  flames.  The  numerical  simulations  were  conducted  using  the  updated  multi-component 
diffusion  model  and  all  four  kinetic  models.  The  results  reported  for  both  (|)  =  0.777  and  0.965 
exhibit  very  similar  behavior  to  those  of  Fig.  15,  in  which  the  updated  multi-component  diffusion 
also  was  used.  For  low  Tad's  (high  N2-dilution)  GRI30  under-predicts  the  experimental  data  by  a 
factor  of  ~2  and  over-predicts  them  at  high  Tad's.  For  both  (|)'s  MKYD99  predicts  closely  the  low- 
Tad  data  and  over-predicts  the  high-Tad  data  by  a  factor  of  1.5-2.  The  predictions  of  TRM04  and 
OPM04  are  consistent  with  those  of  Fig.  15  at  the  corresponding  Tad's. 
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Adiabatic  Flame  Temperature,  K 

1687  1407  1221  1687  1407  1221 


Fig.  17.  Experimental  and  numerical  for 
N2-diluted,  ^  =  0.777  H2/air  flames. 

The  effect  of  reaction  model  is  shown. 


Adiabatic  Flame  Temperature,  K 
1616  1403  1249  1616  1403  1249 


N2/ Air  Molar  Ratio 

Fig.  18  Experimental  and  numerical 
for  N2-diluted,  ([)  =  0.965  H2/air  flames. 
The  effect  of  reaction  model  is  shown. 


Fig.  19.  Logarithmic  sensitivity  coefficients  of  reactions  and  mass  diffusion  coefficients  onK^f 
for  ultra-lean  and  N2-diluted  near-stoichiometric  H2/air  flames  calculated  using  the  TRM04 
model. 

Figure  19  depicts  the  logarithmic  sensitivity  coefficients  of  on  the  kinetics  and  mass 
diffusion  coefficients  for  ultra-lean  H2/air  flames  with  <1)  =  0.289  and  ([)  =  0.314.  While  the 
logarithmic  sensitivity  coefficients  for  kinetics  are  at  most  around  the  value  of  one  (Fig.  19a),  the 
corresponding  ones  for  diffusion  (Fig.  19b)  well  exceed  one,  and,  in  particular,  for  the  H2  values 
of  the  order  of  10  to  nearly  30(!)  are  realized. 
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3.4  Non-Premixed  Ignition  ofH2  and  Gaseous  Hydrocarbons 

A  new  approach  was  developed  for  ignition  studies  in  the  counterflow  configuration.  It 
includes  the  use  of  ultra-lean  H2/air  flames  as  the  ignition  source;  Ref  19  contains  the  pertinent 
methodology  and  results  described  in  this  section.  The  hot  excess  air  dominates  the  composition 
of  the  post-flame  products,  while  minor  amounts  of  H2O  are  present.  The  temperature  of  those 
products  can  be  controlled  by  either  varying  the  H2/air  equivalence  ratio  or  by  adding  small 
amounts  of  CO  in  cases  in  which  higher  ignition  temperatures  and  low  strain  rates  are  required. 
Alternatively,  in  cases  in  which  ignition  temperatures  below  1 100  K  were  required  a  platinum 
screen  was  used  at  the  burner  exit  to  catalytically  support  H2/air  flames  at  such  low  temperatures. 
The  temperature  just  downstream  of  the  platinum  screen  was  measured  by  a  thermocouple  and 
was  used  as  boundary  condition  in  the  numerical  simulations.  The  use  of  the  platinum  screen 
was  required  for  the  ignition  of  H2,  C2H4,  and  C2H2  while  the  use  of  H2/air  flames  at  various 
equivalence  ratios  and  with  various  amounts  of  added  CO  is  required  for  the  ignition  of  saturated 
hydrocarbons.  Premixed  ignition  was  achieved  by  counter-flowing  the  hot  products  against  an 
air-jet  in  which  the  fiiel  was  slowly  added  to  the  air  stream.  Non-premixed  ignition  was  achieved 
by  counter-flowing  the  hot  products  against  a  N2  jet  in  which  the  fuel  was  slowly  added.  The 
non-premixed  ignition  of  H2,  CH4,  C2H6,  C2H4,  C2H2,  and  CsHg  has  been  studied.  Figure  20 
depicts  the  effect  of  H2  and  C2H4  addition  on  the  ignition  of  CH4  for  a  global  strain  rate. 
Kg  =  44  s  '.  These  data  were  determined  by  using  a  (|)  =  0.30  H2/air  flame  as  the  ignition  source, 
and  by  adding  various  amounts  of  H2  and  CH4  in  various  mixtures  of  C2H4/N2.  These  results  are 
of  relevance  to  the  ignition  response  of  the  products  of  thermal  cracking  of  endothermic  fuels 
[e.g.  2].  The  ignition  promotion  caused  by  H2  and  C2H4  is  apparent.  Numerical  studies  showed 
that  no  mechanism  could  predict  the  experimental  results  shown  in  Fig.  20. 


Nonpremixed  CH4  ignition  over  H2  flame  with  the  addition  of 
C2H4  and  H2  additives,  Kp-40  s*^ 
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Fig.  20.  Non-premixed  ignition  of  gaseous  fuels 


Fig.  21.  Premixed  ignition  of  liquid  fuels 
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3.5  Premixed  Ignition  of  Liquid  CH3OH,  C2H5OH,  n-CjHjs,  and  iso-CsHjg  Fuels 

The  premixed  ignition  of  practical  liquid  fuels  has  been  also  studied.  Representative  results 
are  shown  in  Fig.  21  and  have  been  achieved  by  increasing  the  temperature  of  a  (IfO.256  H2/air 
flame  through  CO  addition.  As  expected,  the  higher  was  the  temperature  of  the  H2/CO/air  flame 
the  lower  was  the  ignition  equivalence  ratio  of  the  test  fuel.  The  greater  ignition  propensity  of  the 
alcohol  fuels  as  well  as  that  of  iso-CgHig  compared  to  n-C7Hi6  is  notable. 
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